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Abstract 

We use first-principles calculations to study the effect of magnetic disorder and electron corre- 
lations on the structural and thermodynamic properties of CrN. We illustrate the usability of a 
special quasirandom structure supercell treatment of the magnetic disorder by comparing with co- 
herent potential approximation calculations and with a complementary magnetic sampling method. 
The need of a treatment of electron correlations effects beyond the local density approximation is 
proven by a comparison of LDA+U calculations of structural and electronic properties with exper- 
imental results. When magnetic disorder and strong electron correlations are taken into account 
simultaneously, pressure and temperature induced structural and magnetic transitions in CrN can 
be understood. 
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I. INTRODUCTION 



Transition metal nitrides have attracted much interest due to their excellent performance 
in a long range of industrial application such as hard protective coatings on cutting tools, 
diffusion barriers, and wear resistant electrical contacts. CrN is not as hard as for instance 
TiN [1], but it is superior to TiN in solving large concentrations of AIN in the rock salt 
phase giving rise to the Cri.xAUN solid solutions highly valued in hard coatings applica- 
tions [2-4]. Furthermore, CrN on its own can be found in coating application for metal 
forming and plastic moulding purposes [5, 6]. From a fundamental physics point of view 
the study of CrN has provided new insights but also raised questions about magneto-driven 
structural transitions. It is known that a magnetic order-disorder transition at temperatures 
around 280 K is associated with an orthorhombic to cubic structural transition [7], recently 
observed to be reversible under small pressures [8]. On the other hand in epitaxially stabi- 
lized cubic thin films no sign of magnetic ordering has been seen [9, 10]. On the theoretical 
side the concept of magnetic stress has been introduced and used within a local density 
approximation (LDA) framework to explain the orthorhombic distortion [11, 12]. Recently 
it was shown that taking strong electron correlations into account in the calculations at 
the level of the local spin density approximation plus a Hubbard U-term (LDA+U) could 
improve the agreement between calculations and experiments by opening up a small band 
gap at the Fermi level [13]. However, all these calculations considered only ordered mag- 
netic structures while most experimental measurements, especially of the band structure, 
are performed above the Neel temperature. 

Unlike what is sometimes assumed, most magnetic systems retain magnetic moments 
also above their critical Curie or Neel temperature. Indeed local moments are typically 
present although long-range order between them is lost. CrN is such a system where the ex- 
perimentally observed structural (lattice spacing) and electronic properties (semiconducting 
behavior) of the paramagnetic cubic phase can not be even qualitatively reproduced by non- 
magnetic calculations [11]. At the same time, when performing first-principles calculations 
modelling such disordered cases, ordered magnetic structures should not be used because 
they might give rise to order-specific features, like the well known mixing anomaly in the 
Fei-xCix [14] system. This means that a disordered magnetic state must be considered in 
order to fully understand the physics of paramagnetic CrN at room temperature. 
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For such a purpose the disordered local moments (DLM) [15] method has been suggested 
and implemented within the coherent potential approximation (CPA) [16] treatment for 
disorder. The CPA-DLM method was used in Ref. [17] to demonstrate the importance of 
the magnetic degree of freedom when paramagnetic cubic CrN was alloyed with AIN. Even 
though the DLM-CPA treatment of magnetic disorder is an excellent approximation in many 
cases, a direct method for calculating the electronic structure of magnetically disordered 
systems within a conventional supercell methodology is highly desirable. This is so since the 
CPA is most often combined with other approximations, e.g. the spherical approximation 
for the single-particle potential, imposing certain limitations on the treatment of materials 
with complex underlying crystal lattices. Moreover, in magnetic alloys, such as Cri.j-Alj-N 
or Fei_2;Nia; [18, 19], or in the presence of defects such as nitrogen vacancies in CrN, local 
lattice relaxations and other local environment effects might be important and they are 
beyond the reach of the single-site CPA theory. Furthermore, if the magnetic and vibrational 
thermodynamics of solids are ever to be treated simultaneously on the same footing, a 
supercell treatment of magnetism, compatible with quantum molecular dynamics simulations 
needs to be developed. 

In this work we apply two different supercell approaches to treat the magnetic disorder 
of the paramagnetic phase of CrN and compare them with CPA-DLM calculations. Firstly 
we apply the special quasirandom structure (SQS) method [20], developed to treat chemi- 
cally disordered alloy systems. Secondly, to gain further confidence, we propose a magnetic 
sampling method (MSM) and show that the two methods give equivalent results. 

Furthermore we investigate the impact of strong electron correlations on structural and 
electronic properties of CrN at the level of LDA-I-U calculations. Considering both magnetic 
disorder and strong electron correlations simultaneously we analyze the magneto-structural 
transition in CrN. 

II. CALCULATIONAL DETAILS 

In this work electronic structure calculations are performed within a density functional 
theory (DFT) framework and the projector augmented wave (PAW) method [21] as imple- 
mented in the Vienna ab-initio simulation package (VASP) [22, 23]. Both the local spin 
density approximation (LDA) [24], the generalized gradient approximation (GGA) [25], and 
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a combination of the LDA with a Hubbard Coulomb term (LDA+U) [26, 27] methods are 
used for treating electron exchange-correlation effects. The Hubbard term is applied only 
to the Cr 3d orbitals. In this implementation of the LDA+U method there is only one 
free parameter corresponding to U'^-^''^=(U-J). In the following the simple notation U is used 
for U^-^-^. The energy cut-off for plane waves included in the expansion of wave-functions 
are 400 eV. Sampling of the Brillouin zone was done using a Monkhorst-Pack scheme [28] 
on a grid of 5x5x5 (64-atom supercells), 9x9x9 (48-atoms supercells), 13x13x13 (8-atom 
cells), and 21x21x21 (2-atoms cells) k-points. To find the optimal cell geometry for the or- 
thorhombic structure an automatic optimization procedure was used independently for each 
volume. We also apply the exact muffin-tin orbitals (EMTO) method [29, 30] including the 
full charge-density technique [31] in which the CPA-DLM treatment of magnetic disorder is 
implemented. The EMTO basis set included s, p, d, and f-orbitals and the total energies 
were converged within 0.5 meV/f.u. with respect to the density of the k-point mesh. 

III. SUPERCELL APPROACH TO MAGNETIC DISORDER 
A. The energy of a disordered magnet 

In this work we discuss the thermodynamics of the high temperature paramagnetic state 
of CrN within the classical Heisenberg description, known to work well for many itinerant 
magnetic systems. This state can be described as a disordered distribution of Cr magnetic 
moments on the lattice, lacking long-range order. Similar situations are believed to be 
present in many other systems such as bcc Fe and its alloys at high temperatures [32, 33], 
NiMnSb above the Curie temperature [34] as well as a large number of f-electron systems [35]. 
The first problem we encounter is thus how to calculate the energy of such a disordered mag- 
net using a supercell technique. In order to create an adequate supercell that can be used 
we need to know the characteristics of the random distribution. An ideally random distri- 
bution of magnetic spins, corresponding in the Heisenberg model to infinite temperature, is 
characterised by the vanishing of all the average spin correlation functions 
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where a corresponds to a specific coordination shell, ej is a unit vector in the direction of 
the magnetic moment on site i (e^ = Si/M where M is the magnitude of the magnetic 
moments) and N is a normalisation constant. If the magnetic properties are approximated 
within a classical Heisenberg model hamiltonian, we get 



where is the number of atoms in the a:th coordination shell on the lattice. 

The magnetic moments in CrN have been shown to be formed due to a magnetic split 
of the Cr t2g non-bonding d-states present at the Fermi level [11]. Since large local Cr- 
moments are present and showing rather stable values regardless if they are ordered in a 
ferromagnetic configuration, antiferromagnetic configuration or a disordered local moments 
configuration [17], it seems reasonable that a Heisenberg model hamiltonian could describe 
the magnetism of CrN rather well. We note that this model, although limited to the first two 
nearest neighbour interactions, was applied in Ref. [12] for analyzing magnetic induced stress. 
However, one should be aware that it has been shown that the interaction parameters Jij 
could depend quite substantially on the global magnetic state even if the magnetic moments 
are almost constant [34] , underlining the importance of a reliable method directly accessing 
the disordered magnetic state. 

Two important observations can be made from Eq. 2. Firstly, in order to calculate the 
total energy of the ideal random distribution, the structure we use in simulations has to fulfill 
the condition in Eq. 1 only for the coordination shells a where the interaction parameters 
Ja are non- negligible. Secondly, one realizes that even though the disordered spin state 
is utterly non-coUinear, its energy can be calculated using a collinear state as long as the 
parallel (Sj • Sj = +1 • M^) and anti-parallel (Sj • Sj = —1 • M^) spin pairs exactly cancel on 
each relevant coordination shell. Thus, it is these characteristics that we should aim for in 
our simulation of the high temperature paramagnetic state. 

It is not directly obvious how to create a supercell fulfilling these properties but it is in 
fact a situation very similar to the problem of modeling chemical disorder in the form of 
a binary substitutional random alloy [36]. The DLM-CPA method mentioned above can 
actually be seen as just a CPA treatment of a random alloy of atoms with spin up and 
spin down oriented magnetic moments. If we follow this analogy to the supercell framework 
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it is the special quasirandom structure (SQS) methodology, first suggested by Zunger et 
al. [20], that has proven to be the most accurate approach for direct calculations of the total 
energy or related properties of the disordered state. The agreement between the CPA and 
SQS methods for chemical disorder in the Bl structure was demonstrated for the Tii_^.Al^.N 
system in Ref. [37] using the same electronic structure methods as in this work. Here we use 
the 48-atom (24 metal + 24 nitrogen) SQS structure suggested in Ref. [37] for Tio.5Alo.5N 

't I 

to model the CrQ^Crp^N paramagnetic phase. That structure has < $q, >= for the 
first 7 coordination shells with the exception of a small non-zero value on the fifth shell. 
For comparison reasons we also use a 64-atom (32 chromium + 32 nitrogen) SQS-structure 
based on a cubic 2x2x2 conventional unit cell geometry with < $c» >= for the first 6 shells 
with the exception of a small non-zero value on the third shell. We used the SQS method 
in a study of the bulk modulus of the paramagnetic phase of CrN [38] as well as to get a 
CrN reference energy in a study of Tii^^^Crj-N [39]. However the reliability of the method 
was questioned in Ref. [40]. 

To check if the SQS method is reliable to model magnetic disorder on the cubic lattice of 
CrN we show in Fig. 1 a comparison of energy-lattice parameter curves for different mag- 
netic states including the disordered state calculated with the SQS-PAW (left panel) and 
CPA-EMTO (right panel) methods employing the GGA functional for exchange- correlation 
effects. The energies are given relative to the non-magnetic energy minimum. The agree- 
ment between the two different treatments of disordered magnetism is clearly seen as the 
disordered state is placed in a very similar relation to ordered magnetic and non-magnetic 
calculations in the two methodological frameworks. 

Here we note that the accuracy of the DLM-CPA treatment of a completely disordered 
magnetic state is established analytically on the single-site level in Ref. [15]. We therefore 
view very good agreement between SQS and DLM-CPA calculations as a strong proof that 
the former technique is capable to describe the energy of a paramagnetic state, at least at 
the same level of accuracy. In particular, it is clear that the SQS method does not suffer 
from imposed periodic boundary conditions or the fact that one technically speaking deals 
with one selected antiferromagnetic configuration. As soon as conditions given by Eqs. 1 
and 2 are fulfilled, the SQS represents a quasi-random rather than ordered magnetic state. 
As for some minor differences between the results presented in the two panels of Fig. 1, they 
come from the usage of different underlying methods for the electronic structure calculations. 
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FIG. 1: (Color online) The energies as a function of lattice parameter of different magnetic states of 
cubic Bl CrN as calculated with the PAW and EMTO methods and the GGA exchange-correaltion 
functional. The disordered phase (circles) is modelled with the SQS method in the PAW calcula- 
tions and the DLM-CPA method in the EMTO calculations. 

PAW and EMTO. 

In order to further establish the reliability of the SQS approach, using the same PAW 
methodology, we first compare the results calculated for the two different SQS geometries 
considered in this work. The energy difference between them is 0.003 eV/f.u. Due to the 
translational symmetry the SQS based on the 2x2x2 conventional unit-cells has the problem 
that the correlation function on the 8:th correlation shell is exactly 1 and this is probably 
the main source of the small difference between the two SQSs. 

Next we suggest a different method to calculate the energy of a magnetic state approxi- 
mating that in Eq. 1.: The magnetic sampling method (MSM). Within the MSM the direc- 
tions (up and down) of magnetic moments of the Cr-sites of a large supercell are chosen using 
a random number generator. A large set of different such distributions, magnetic samples, 
are then created. Their energies are calculated and the average energy is taken as the energy 
of the disordered state. Individually these supercells typically do not satisfy the conditions 
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FIG. 2: (Color online) The calculated energies of the different random number generated magnetic 
configurations of a 32 Cr-atoms Bl CrN supercell (circles). The accumulated average of the random 
number generated configurations, called magnetic iterations, is shown with a solid line. All values 
are plotted relative the energy minimum of a SQS based magnetic configuration on the same 
underlying structure. 



of Eqs. 1 and 2 but their average should, given that a sufficiently large number is consid- 
ered. Although the SQS formalism was suggested in a reaction against the inaccuracies of 
random number distribution schemes [20], we note that the supercell sizes and particularly 
the number of calculations possible to treat with todays computational resources are orders 
of magnitude larger as compared to those back in 1990. 

Fig. 2 shows the calculated energies of 40 different randomly generated magnetic samples 
with up and down collinear moments on the ideal lattice points of a 32 Cr-atoms (2x2x2 
conventional unit cells) Bl CrN supercell. Their accumulated average energy is shown with 
a solid line and compared to the energy of the SQS generated configuration on the same 
underlying geometry. The latter is taken as the reference energy. Although the individual 
energies of the randomly generated supercells can differ as much as -0.02 eV/f.u. and -1-0.035 
eV/f.u. from the SQS-value, the accumulated average after the consideration of 40 different 
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samples, called magnetic iterations, differs only by 0.001 eV/f.u. Between iteration 20 and 
40, the average is never more than 0.001 eV/f.u. above or 0.0015 eV/f.u. below the value 
at iteration 40 showing that the mean energy is converging. If needed one could use even 
more sampled structures to converge the value with a higher accuracy than the 40 used 
in the present work. With this said we conclude, by our comparison with the CPA-DLM 
calculations as well as the internal agreement between SQS and MSM methods that both 
the considered supercell approaches can be used to calculate the total energy of a disordered 
collinear magnetic state of CrN with an accuracy of a few meV/f.u. on a fixed ideal Bl 
lattice. 

B. Non-collinear considerations 

Finally, we note that for systems where the energetics of the non-collinear disordered state 
is believed not to be well described by Eq. 2, for instance due to non- negligible contributions 
from bi-quadratic terms in the Hamiltonian: 

Hmag = - ^ Jij^i " " ^ Kij ■ Bjf = - ^ ( Ja^a < $q > +Kana < ^„ >) , (3) 

where in the fully disordered state 

iJCct 

the MSM method could still be used with a straight forward generalisation: Instead of ran- 
dom number generation of collinear spins (up and down in z), one can generate a set of dif- 
ferent non-collinear supercells with six types of local moments describing up and down along 
X, y, and z. Such a set, given that the supercells are large enough and that the number of 
supercells is large enough, will reproduce both the bi-linear (Eq. 1), and bi-quadratic (Eq. 4) 
correlation functions of the disordered state. In principle, but cumbersome in practise, also 
the SQS method can be used in this case by constructing a large supercell of 6 compo- 
nents with vanishing correlation functions between them all. Unfortunately, first-principles 
calculations of non-collinear magnetic systems are rather time consuming, nevertheless we 
have performed non-collinear MSM calculations within the LDA+U approximation for the 
volume obtained with the collinear approximation and will get back to this issue below. 
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IV. EFFECT OF STRONG ELECTRON CORRELATIONS 



Having established methods to treat the magnetic disorder within a supercell framework 
we now turn to the problem of electron exchange-correlation energies. Even though LDA 
calculations qualitatively revealed the energetics of CrN [11, 12] the electronic structure 
did not reproduce the experimentally observed semiconducting behavior. Thus one could 
doubt the accuracy of LDA predictions of structural and magnetic energy differences of 
relevance for understanding the orthorhombic to cubic transition in this system. Recently 
Herwadkar et al. studied CrN using the LDA+U approach with focus on the electronic 
structure of ordered magnetic structures. They calculated the value of U and J, the screened 
Coulomb and exchange terms respectively, using constrained LDA approaches and achieved 
U=3 eV and J=0.9 eV [13]. However, they suggested a span of U values from 3-5 eV to 
be reasonable. Even though such an ah-inito approach to obtain the values of the U and J 
parameters are appealing, the uncertainties are to large for a quantitative a thermodynamic 
analysis. 

In order to obtain the most suitable value U^^^ , for which the LDA+U method best 
describes the relevant properties of CrN, we perform a careful comparison of structural 
parameters and electronic structure obtained with LDA+U calculations for various values of 
U with experiments. For comparison also the results obtained with the generalized gradient 
approximation, GGA, are presented. 

First the lattice parameter of the cubic paramagnetic phase, modeled with the SQS ap- 
proach, is presented in Fig. 3. The experimental value obtained for bulk CrN by Corliss et 
al. is 4.13 A [7], while both Herle et al. [41] and, more recently, Rivadulla et al. [8] obtained 
4.148 A. Values obtained for CrN in thin films are typically slightly larger [9, 10] but include 
strain effects not considered in the calculations. The result obtained with the pure LDA 
functional is a=4.022 A. Using GGA we obtain a=4.149 A. The LDA+U approach gives 
increasing lattice spacing with increasing U. The lattice parameter of Corliss et al. is ob- 
tained with U=2.9 eV. The value of Herle et al. and Rivadulla et al. is obtained with U=3.8 
eV. Since the reported experimental values are measured at room temperature while the 
calculations, with the exception for the magnetic disorder, corresponds to a K situation, 
one might object that the comparison is not completely fair. However, thermal expansion 
between K and room temperature is typically small in this class of hard ceramics. 



10 




1 2 3 4 5 6 



U(eV) 

FIG. 3: (Color online) The calculated cubic lattice parameter for the magnitically disordered state 
within the LDA+U (circles) and GGA (square) approximations. The experimental bulk value 
found by Corliss et al. [7] and the common value found by Herle et al. [41] and Rivadulla et al. [8], 
are shown with solid and dashed horizontal lines respectively. 

Since also the physics of the orthorhombic phase must be well described by our theoretical 
model, we compare in Fig. 4 the calculated value of the angle a with the experimental finding 
in Refs. [7, 8]. a describes the angle between the axis of the conventional unit cell of the 
cubic Bl lattice, that is distorted in the orthorhombic phase, see the inset in Fig. 4. The 
experimental value is 88.3°according to Corliss et al. [7] or 88.4° according to Rivadulla et 
al. [8], while both LDA and GGA calculations underestimate this angle, thus overestimating 
the distortion. On the other hand, as can be seen in Fig. 4, a value in agreement with the 
experiment is obtained with the LDA+U method for U=3.0 eV. 

Finally, in Fig. 5 we compare the calculated total electronic density of states of the valence 
band of the disordered magnetic state (calculated with the SQS method) with the ultravio- 
let photoemission spectroscopy measurement of the cubic paramagnetic phase obtained by 
Gall et al. [9]. Figure 5 shows in different panels (from top to bottom) the DOS obtained 
with GGA, and LDA+U with U from eV (LDA) to U=5 eV. In all panels the experimental 
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FIG. 4: (Color online) The calculated angle a in the distorted orthorhombic state between the 
axis of the conventional Bl-cell. Values using the LDA+U approximation (circles) and the GGA 
approximation (square) are shown together with the experimental value found by Corliss et al. [7] 
(solid horizontal line) and Rivadulla et al. [8] (dashed horizontal line). Inset: The definition of the 
angle a with respect to the orthorhombic unit cell (dashed lines) shown in a 001-plane of the cubic 
Bl structure. The relative directions of the spins of Cr atoms following Ref. [7] are also shown 
(solid circles: z=0.0, open circles: z=0.5). 

results are shown with dashed hnes. In all cases we use the lattice spacing corresponding 
to the equilibrium of the particular choice of exchange-correlation scheme. The GGA, and 
even more so the LDA, gives an overlap of the peaks close to the Fermi level. These peaks 
correspond primary to Cr spin-up non-bonding (below Ep) and a combination of spin down 
Cr non-bonding and Cr spin up anti-bonding states (above Ep). Both peaks also have small 
ad-mixture of N p-character [13]. When the U value in the LDA-I-U approach is increased, 
the occupied Cr spin up non-bonding state becomes more localized and shift down in energy. 
At the same time the unoccupied states are shifted up and at a value of U between 2 and 3 
eV a small gap opens at Ep, in agreement with the experiment [9]. Actually the LDA-I-U 
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FIG. 5: (Color online) The calculated valence band electronic density of states (solid line in all 

frames) of the cubic magnetically disordered state using the GGA approximation (top panel) 

and LDA+U approximation with different values of U. For comparison the experimental UPS 

spectroscopy measurement by Gall et al. [9] is shown by a dashed line in all frames. 
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approximation with U=3 eV excellently describes the Cr spin up non-bonding state, while 
the bonding states on the other hand are shifted to slightly too high energies. 

In summary, the LDA+U approximation with U- values between about 3 and 4 eV repro- 
duces the cubic paramagnetic lattice parameter, the value U=3 eV reproduces experimental 
measurement of the angle a in the orthorhombic phase, and a value of U between 2 and 
3 eV gives a good description of the electronic structure of the valence band of the cubic 
paramagnetic phase. Thus we conclude that within the LDA+U approximation, the value 
U^-^^=3 eV gives an optimal description of the physical properties of the system. This value 
is safely within the range suggested in Ref. [13] corresponding to U=3.9 if their choice of 
J=0.9 is used. 

In the following sections we use the LDA+U approximation with effective U=3 eV in our 
calculations. However, since the GGA approximation has already been used in many studies 
on related systems, we present also results using the GGA to illustrate the effect of strong 
electron correlations. 

V. ENERGETICS OF MAGNETIC AND CRYSTALLOGRAPHIC PHASES OF 
CrN 

We now have the theoretical tools needed to study the magneto-structural transition in 
CrN. In Fig. 6 we consider the total energies, as a function of volume of different crystal- 
lographic and magnetic phases of CrN. The left panel shows the results obtained with the 
GGA functional while the right panel shows the results obtained with LDA+U (U=3 eV). 
Cubic phases are shown with open symbols while orthorhombic (orth.) structures are shown 
with solid symbols. The paramagnetic cubic phase is modeled with the SQS disordered 
local moments method (denoted cubic dim) keeping the atoms fixed at Bl lattice points. 
The energy of a disordered magnetic configuration on the lattice points of the orthorhombic 
structure is also shown for comparison and denoted orth. dim. The minimum energy of the 
cubic collinear disordered magnetic phase is taken as the reference value. Furthermore, the 
energy of the experimental ground state structure [7], the orthorhombic distorted double 
[Oil] -layered antiferromagnetic ([Oil] 2 afm) structure schematically shown in the inset of 
Fig. 4, is shown as is the energy of the same magnetic ordering on the cubic lattice. Also 
a single layer [001]-ordered antiferromagnetic state ([001]i afm) on the cubic lattice, the 
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FIG. 6: (Color online) The calculated energy versus volume curves for different magnetic states in 
CrN using the GGA approximation (left panel) and LDA+U (U=3eV) (right panel). The energy 
minimum of the cubic collinear disordered local moment state (DLM) is used as reference energy. 

ferromagnetic cubic phase, and in the case of LDA+U the energy of the non-colhnear cubic 
magnetically disordered phase are shown. 

One can see that the GGA and LDA+U calculations are in reasonable qualitative agree- 
ment with each other. However, the GGA gives considerably larger energy differences be- 
tween the orthorhombic and cubic phases as compared to the LDA+U calculations. Also the 
order of the two considered antiferromagnetic states for the cubic phase are reversed. In the 
GGA framework the [001] i antiferromagnetic state is lower in energy, in line with previous 
works [11, 17]. In the LDA+U framework on the other hand, the [011]2 antiferromagnetic 
state is lowest in energy also in the cubic phase. 

One interesting comparison can be made between the introduction of a Hubbard U- 
term in this work and the alloying of GrN with AIN studied in Ref. [17]. In that work 
it was found that the DLM-CPA state became lower in energy as compared to the [001]i 
antiferromagnetic state when a certain amount of Al was substituted for Gr. The point is 
that upon alloying of GrN with AIN [17], and more generally upon alloying transition metal 
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nitrides with AIN [42], the inclusion of Al favours a locahsation of the transition metal non- 
bonding d-states. Obviously, the strong electron correlations leads to a similar effect which 
explains the similar evulotion of the magnetic energies in the two cases. 

VI. NON-COLLINEAR MAGNETISM IN CrN 

Using the LDA+U approximation the impact of non-collinear magnetism on the energy 
of the disordered cubic phase has been investigated with the MSM method by averaging 
over 35 different random-number generated non-collinear distributions, a procedure capable 
of reproducing the condition of both Eqs. 1 and 4 for the relevant coordination shells. The 
average energy is found to converge in a similar manner as the coUinear case described above. 
The energy of the non-collinear disordered magnetic state is found to be 0.010 eV/f.u. below 
that of the collinear disordered state and is marked in Fig 6 with a striped circle. This energy 
difference is not particularly large with respect to e.g. the orthorhombic antiferromagnetic 
state or the cubic ferromagnetic state, but as will be seen below, do have a quantitative 
effect on the simulation of the cubic-to-orthorhombic transition temperature. In our case 
the improved level of theoretical modelling when including non-collinear effects, came to the 
cost of an approximately seven-fold increase in computational time with respect to otherwise 
identical settings for the calculations of collinear disordered structures. 

VII. THE MAGNETO-STRUCTURAL TRANSITION IN CrN 

Qualitatively our calculated values agree with the experimental observation of the stabil- 
ity of the orthorhombic state at low temperatures and a cubic state at higher temperatures. 
In this work we denote the temperature for this structural transition Ts in order not to 
confuse it with the hypothetical isostructural Neel temperatures, Tjv, of a magnetic order- 
disorder transiton on a fixed lattice. This is so since the disordered paramagnetic phase has 
a considerable magnetic entropy making it more competitive at higher temperatures. Since 
the energy of the cubic dim-state is considerably lower than the orthorhombic dim-state the 
magnetic disordering is accompanied by a structural transition. The fact that no signs of 
magnetic ordering was observed in the epitaxially stabilized cubic phases in Ref. [9, 10] can 
be understood from the fact the the energy differences between antiferromagnetic and disor- 
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dered magnetic phases in the cubic geometry are small, especially if non-coUinearity is taken 
into account, indicating a very low Neel temperature for undistorted cubic CrN. In the or- 
thorhombic structure on the other hand the difference is almost an order of magnitude larger 
indicating that the orthorhombic antiferromagnetic state should be well below its isostruc- 
tural Neel temperature at the experimental transition point T^^^* = 280 — 287 K [7, 8]. 
This result is in line with the experimental observation that the extent of the orthorhom- 
bic distortion, measured with the value of the angle a, is in principle the same at 286 K: 
88.4 ° [8], 273 K: 88.3 ° [43], and 77 K: 88.3 ° [7] where the authors of the latter reference 
stated that no differences was seen when the temperature where further decreased down to 
the liquid Helium regime. If there had been a large degree of partial magnetic disorder in 
the orthorhombic phase one would expect a change in the value of this angle. Thus it is 
reasonable to assume that the transition is governed by the competition in terms of free 
energy between a disordered paramagnetic cubic phase with high magnetic entropy and a 
highly ordered antiferromagnetic orthorhombic phase with low magnetic entropy. Such a 
phase transition, including an abrupt change in both energy and entropy is in line with the 
experimental finding of a first order phase transition displaying a hysteresis behavior during 
heating-cooling cycles [43]. 

Using our obtained structural energy differences we can estimate the transition tempera- 
ture theoretically. At temperatures considerably above the (isostructural) Neel temperature, 
such as in the case of paramagnetic CrN at room temperature, the entropy of a Heisenberg 
system is well described by the mean-field term 

3""^ = kBln{M+l), (5) 

where M is the magnitude of the magnetic moment (in units of fis) and ks is the Boltz- 
mann constant. In the LDA+U approach we find that the average magnetic moments are 
]^LDA+u ^ 2.82/iB- In the GGA calculation M^^^ = 2.49/iB. 

Using these values and the approximations above the transition temperature can be 
obtained from the condition that at the critical temperature, the two phases should have 
the same free energy, F: 

F'^JlHiTs) = F;:1{Ts) ^ EZft = EZta - TsS^^ (6) 
where Ts denotes the critical temperature for the structural transition which is in re- 
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ality also the magnetic ordering temperature. In our case we get with the colhnear ap- 
proximation for mag netic disorder ^"irpLDA+u ^ ^gg ^ ^-^j^ 

non-coUinear calculations 

ncrpLDA+u ^ ^ ^j^j^ ^Yie collinear GGA calculations = 1030 K. This mean- 

field estimates should be compared to the experimental value of = 280 — 287 K [7, 8] . 
It is well known that the mean field approximation could overestimate magnetic ordering 
temperatures by as much as 50% as compared to more reliable thermodynamics treatments. 
However, in the GGA calculation the error is so large that we instead interpret this result as 
one more argument that the GGA approximation overemphasize the orthorhombic distor- 
tion, both in geometric distortions visible in Fig. 4, and in the energy differences between the 
orthorhombic and cubic phases. The transition temperature derived from the non-coUinear 
LDA+U calculation is closer but still considerably above the experimental measurement, 
giving an overestimation of Ts with 44%. On the other hand, the large relative error may be 
due to the low absolute value of the transition temperature. The absolute error in Ts is just 
a little bit above 100 K for LDA+U treatment including non-coUinear effects, which is the 
most typical case for first-principles simulations of structural phase transitions. In particular 
we suggest that magnetic short-range order coupled with the vibrational degree of freedom 
is of importance to quantitatively determine the transition temperature in CrN. It should 
be noted that also nitrogen off-stoichiometry could influence the transition temperature. 
Moreover, one should remember that the LDA+U approach is an approximate method not 
free from errors, for instance in the exact choice of U, possibly affecting quantitative values 
of the important structural energy difference. 

Recently RivaduUa et al. showed that the temperature induced orthorhombic to cubic 
phase transition could be reversed with increasing pressure [8]. At room temperature, a 
pressure as low as 1-2 GPa was enough to push the system back into the orthorhombic 
structure. Magnetic measurements showed that it was the antiferromagnetic ordered struc- 
ture that reappeared [8] . Qualitatively the pressure effect can be understood from the results 
in Fig. 6: Since the orthorhombic phase is slightly lower in volume as compared to the para- 
magnetic cubic phase, it will be relatively more favorable at elevated pressures according to 
the minimization of the Gibb's Free energy 



G{T,p) = E + pV -TS. (7) 
The derived pressure-temperature phase diagram of CrN is shown in Fig. 7. The results 
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FIG. 7: (Color online) The calculated pressure-temperature phase diagram of CrN using the 
LDA+U (U=3eV) method and considering the disordered magnetic phase as collinear (green dashed 
line) or non-collinear (green diamond). The experimental measurements of Rivadulla et al. [8] at 
low pressures are shown with a bold black line while a linear extrapolation to higher pressures is 
shown with a thin black line. 

from the theoretical calculations using Eq. 7 with the mean field approximation for the 
magnetic entropy of the cubic phase, Eq. 5, and the LDA+U approximation for exchange- 
correlation effects are compared to the experimental low-pressure results from Ref. [8] and 
a linear extrapolation of these values to medium pressures. The qualitative picture, with 
increasing transition temperature with increasing pressure is rather well reproduced by the 
calculations although the absolute values of the temperatures are to high. Partially they are 
corrected by the consideration of non-collinear effects. 

These results inspire us to propose a possible way to conclude the discussion of the value 
of the bulk modulus of cubic CrN [8, 38, 40]. We suggest that in order to measure the 
compressibility of the cubic paramagnetic phase of CrN with higher accuracy as compared 
to Ref. [8], the experiment should be conducted at slightly higher temperatures where the 
cubic phase is stable over a larger pressure range. 
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VIII. CONCLUSIONS 



We have used two different supercell approaches to model disordered magnetism of para- 
magnetic materials, the special quasirandom structure method and the magnetic sampling 
method, and applied them for the study of CrN. The SQS method is valid and swiftly ap- 
plicable when non-collinear effects can be neglected. The MSM and SQS mehods are shown 
to give equivalent results for coUinear calculations of disordered local moments on a fixed 
Bl lattice in cubic CrN and both of them agree with DLM-CPA calculations. We show that 
it is straight-forward to extend the MSM method to calculations of non-collinear disordered 
magnetism. The latter is shown to have a quantitative influence on the calculated structure 
transition temperature in CrN. 

CrN is a correlated material which is better described with a LDA+U aproach then 
with the GGA or LDA functionals. By comparing the calculated structural parameters and 
electronic structure of CrN with experiments we find that U'^^^=3 eV is a suitable value to 
use in the simulations. 

Considering both magnetic disorder effects and strong electron correlations, the or- 
thorhombic to cubic phase transition of CrN as a function of temperature and pressure 
can be qualitatively explained. In particular we show that it should be understood as a 
transition from a magnetically ordered orthorhombic phase to a magnetically disordered 
cubic phase. Considering magnetic entropy within the mean field approximation, the cal- 
culated transition temperature T5 = 413 K is a. slight overestimation of the experimental 
value 280 — 287 K. Magnetic short range order coupled with vibrational effects are likely 
to be of importance for determining the quantitative value of T5. Since the transition also 
depends sensitively on the structural energy difference between the cubic and orthorhombic 
phases, which is shown to be very sensitive to the exchange- correlation functional, a strong 
electron correlations method beyond the LDA+U approach might be needed to reveal the 
details of the CrN phase transition. 
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